Extremal transmission at the Dirac point of a photonic band structure 
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We calculate the effect of a Dirac point (a conical singularity in the band structure) on the trans- 
mission of monochromatic radiation through a photonic crystal. The transmission as a function 
of frequency has an extremum near the Dirac point, depending on the transparencies of the in- 
terfaces with free space. The extremal transmission To — ToW/L is inversely proportional to the 
longitudinal dimension L of the crystal (for L larger than the lattice constant and smaller than the 
transverse dimension W). The interface transparencies affect the proportionality constant To, and 
they determine whether the extremum is a minimum or a maximum, but they do not affect the 
"pseudo-diffusive" 1/L dependence of To. 
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I. INTRODUCTION 



In a two-dimensional photonic crystal with inversion 
symmetry the band gap may become vanishingly small 
at corners of the Brillouin zone, where two bands touch 
as a pair of cones. Such a conical singularity is also re- 
ferred to as a Dirac point, because the two-dimensional 
Dirac equation has the same conical dispersion relation. 
In a seminal work Raghu and Haldane investigated 
the effects of broken inversion symmetry and broken time 
reversal symmetry on the Dirac point of an infinite pho- 
tonic crystal. Here we consider the transmission of ra- 
diation through an ideal but finite crystal, embedded in 
free space. 

As we will show, the proximity to the Dirac point is as- 
sociated with an unusual scaling of the transmitted pho- 
ton current / with the length L of the photonic crys- 
tal. We assume that L is large compared to the lattice 
constant a but small compared to the transverse dimen- 
sion W of the crystal. For a true band gap, / would be 
suppressed exponentially with increasing L when the fre- 
quency uj lies in the gap. Instead, we find that near the 
Dirac point / oc 1/L. The 1/L-scaling is reminiscent of 
diffusion through a disordered medium, but here it ap- 
pears in the absence of any disorder inside the photonic 
crystal. 

Such "pseudo-diffusive" scaling was discovered in Refs. 
0, [3] for electrical conduction through graphene (a two- 
dimensional carbon lattice with a Dirac point in the spec- 
trum) . Both the electronic and optical problems are gov- 
erned by the same Dirac equation inside the medium, but 
the coupling to the outside space is different. In the elec- 
tronic problem, the coupling can become nearly ideal for 
electrical contacts made out of heavily doped graphene 
0) 9 ' by suitably matching the Fermi energy in metal- 
lic contacts 0, Q . An analogous freedom does not exist 
in the optical case. 

The major part of our analysis is therefore devoted to 
the question how nonideal interfaces affect the depen- 



dence of / on uj and L. Our conclusion is that 



I/Io - ToW/L 



(1.1) 



at the Dirac point, with Iq the incident current per mode 
and Fo an effective interface transparency. The proper- 
ties of the interfaces determine the proportionality con- 
stant Fo, and they also determine whether / as a function 
of w has a minimum or a maximum near the Dirac point, 
but they leave the 1 /L-scaling unaffected. 

In Sec. |TT] we formulate the wave equations inside and 
outside the medium. The Helmholtz equation in free 
space is matched to the Dirac equation inside the pho- 
tonic crystal by means of an interface matrix in Sec. IIIII 
This matrix could be calculated numerically, for a spe- 
cific model for the termination of the crystal, but to arrive 
at general results we work with the general form of the 
interface matrix (constrained by the requirement of cur- 
rent conservation). The mode dependent transmission 
probability through the crystal is derived in Sec. IIVI It 
depends on a pair of interface parameters for each of the 
two interfaces. In Sec.|V]we then show that the extremal 
transmission near the Dirac point scales oc 1/L regard- 
less of the values of these parameters. We conclude in 
Sec. IVII with suggestions for experiments. 



II. WAVE EQUATIONS 

We consider a two-dimensional photonic crystal con- 
sisting of a triangular or honeycomb lattice in the x-y 
plane formed by cylindrical air-filled holes along the z- 
axis in a dielectric medium (see Fig. [1]). The crystal has 
a width W along the j/-direction and a length L along the 
x-direction, both dimensions being large compared to the 
lattice constant a. Monochromatic radiation (frequency 
uj) is incident on the plane a; = 0, with the electric field 
£'(a;,y)e"^* polarized along the z-axis. 

In the free space outside of the photonic crystal {x < 
and X > L) the Maxwell equations reduce to the 
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FIG. 1: (Color online.) Photonic crystal formed by a di- 
electric medium perforated by parallel cylindrical holes on a 
triangular lattice (upper panel: front view; lower panel: top 
view) . The dashed lines indicate the radiation incident on the 
armchair edge of the crystal, with the electric field polarized 
in the z-direction. 



Helmholtz equation 



(92 + 92) E{x,y) + —E{x,y) = {). 



(2.1) 



The mean (time averaged) photon number flux in the 
x-direction is given by [6| 



Aiftjjj'^ V dx dx 



(2.2) 



Inside the photonic crystal (0 < a: < L) the Maxwell 
equations reduce to the Dirac equation [l| 











idy 



= {uj - lod) 



^2 
^2 



(2.3) 



for the amphtudes ^^i, ^2 of a doublet of two degenerate 
Bloch states at one of the corners of the hexagonal first 
Brillouin zone. 

As explained by Raghu and Haldane [H, 0| , the modes 
at the six zone corners Kp^K'^ (p = 1,2,3), which are 
degenerate for a homogeneous dielectric, are split by the 
periodic dielectric modulation into a pair of doublets at 
frequency lud and a pair of singlets at a different fre- 
quency. The first doublet and singlet have wave vectors 
at the first set of equivalent corners Kp, while the second 
doublet and singlet are at K'p. Each doublet mixes and 
splits linearly forming a Dirac point as the wave vector 
is shifted by Sk from a zone corner. The Dirac equa- 
tion (|2.3p gives the envelope field cx e"^'" '' of one of these 
doublets. 




FIG. 2: Right panels: Hexagonal first Brillouin zone of the 
photonic crystal (top) and dispersion relation of the doublet 
near one of the zone corners (bottom). Filled and open dots 
distinguish the two sets of equivalent zone corners, centered 
at Kp and K'p, respectively. The small circles centered at 
the zone corners are the equal-frequency contours at a fre- 
quency Lu just above the frequency ud of the Dirac point. 
Left panels: Equal-frequency contour in free space (top) and 
corresponding dispersion relation (bottom). A plane wave in 
free space with kx close to fco (arrows in the upper left panel) 
excites Bloch waves in the photonic crystal with k close to 
Ki and K2 (arrows in the upper right panel), as dictated by 
conservation of ky and uj (dotted horizontal lines). 



The frequency lod and velocity vd in the Dirac equa- 
tion depend on the strength of the periodic dielectric 
modulation, tending to ujd = c'\Kp\ = c'\Kp\ — Attc' /3a 
and V]j — c'/2 in the limit of weak modulation. (The 
speed of light c' in the homogeneous dielectric is smaller 
than the free space value c.) 

Eq. ()2.3p may be written more compactly as 



«ud(V • cr)^' = Slu'^, Suj = lo — uj]j, 



(2.4) 



in terms of the spinor ^E* — (5*1, ^'2) and the vector of 
Pauli matrices cr = {(Jx,'^y)- In the same notation, the 
velocity operator for the Dirac equation is ud<t. The 
mean photon number fiux in the .T-direction is there- 
fore given by 



3D = WD**cr.^* = Wr,(*i*2 + *2*l)- 



(2.5) 



The termination of the photonic crystal in the y- 
direction introduces boundary conditions at the edges 
y — and y = W which depend on the details of the 
edges, for example on edges being of zigzag, armchair, 
or other type. For a wide and short crystal, W ^ L, 
these details become irrelevant and we may use periodic 
boundary conditions [^'(a;,0) = ^'(a;, W)] for simplicity. 
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III. WAVE MATCHING 

The excitation of modes near a Dirac point has been 
discussed by Notomi Q, in terms of a figure similar to 
Fig. [21 Because the j/-component of the wave vector is 
conserved across the boundary at x — 0, the doublet 
near Ki = {Kx,Ky) or K2 = {—Kx,Ky) can only be 
excited if the incident radiation has a wave vector k ~ 
(kxyky) with ky near Ky. The conservation of ky holds 
up to translation by a reciprocal lattice vector. We will 
consider here the case of |fe| < \Kp\, where no coupling to 
is allowed. The actual radius of the equal frequency 
contour in the free space at uj = lud will depend on a 
particular photonic crystal realization. 

The incident plane waves -^incident = E^e^^ '^ in free 
space that excite Bloch waves at a frequency 5uj = uj~ujd 
have ky = Ky[l-\-0{5ijj /ujd)] and k^ = ko[l-\-0{5u) /lob)] 
with 

fco = ^{ojdIcY - Kl. (3.1) 

For SliJ <^ uJu we may therefore write the incident wave 
in the form 

i^incidcntCx, y) = E+{x, y)e'''o-+'^yy, (3.2) 

with a slowly varying function. Similarly, the re- 
flected wave will have ky ~ Ky and kx ~ —fco, so that we 
may write it as 

i^rcficctcd(x,y) = ^;_(x,^/)e-^"°^+*^''^ (3.3) 

with slowly varying. 

The orientation of the Brillouin zone shown in Fig. [2] 
corresponds to an armchair edge of the triangular lattice 
at a; = 0. For this orientation only one of the two in- 
equivalent doublets is excited for a given ky. (The other 
doublet at K[, K2 is excited for —ky.) A 90° rotation 
of the Brillouin zone would correspond to a zigzag edge. 
Then a linear combination of the two inequivalent dou- 
blets is excited near ky = 0. For simplicity, we will re- 
strict ourselves here to the case shown in the flgurc of 
separately excitable doublets. 

While the conservation of the wave vector component 
parallel to the boundary determines which modes in the 
photonic crystal are excited, it does not determine with 
what strength. For that purpose we need to match the 
solutions of the Helmholtz and Dirac equations at x = 
0. The matching should preserve the flux through the 
boundary, so it is convenient to write the flux in the 
same form at both sides of the boundary. 

The photon number flux (|2.2[) for the Helmholtz equa- 
tion may be written in the same form as the flux p.5|) 
for the Dirac equation, by 

jn = vh£*(Tx£, (3.4a) 
gpc^fco (. (E++E\ 



(In the prefactor fcp we have neglected corrections of order 
du/iOD.) Flux conservation then requires 

vh£*(Tx£ = VD'i'*crx'i', at a; = 0. (3.5) 

The matching condition has the general form [§] 

* = (wH/wD)^/^Mf, at x^O. (3.6) 

The flux conservation condition p.Sp implies that the 
transfer matrix M should satisfy a generalized unitarity 
condition, 

M-i = axM^ax. (3.7) 

Eq. (|3.7p restricts M to a three-parameter form 

M = e^'^-e^'^^e*"'"- (3.8) 

(ignoring an irrelevant scalar phase factor). The real pa- 
rameters a, /?, 7 depend on details of the boundary at 
the scale of the lattice constant — they can not be deter- 
mined from the Helmholtz or Dirac equations (the latter 
only holds on length scales 3> a). 

We now show that the value of a becomes irrelevant 
close to the Dirac point. At the boundary the incident 
and reflected waves have the form 

•fincidcnt = Eq , iSroficctcd = tEq ' ('^■^) 

with r the reflection coefficient, and Eq = E^{0,y) a 
slowly varying function. Both "spinors" are eigenvectors 
of ax, hence the action of e*"'^= on £ is simply a phase 
factor: 

^incident — ^incident, 

Mfrefiectcd = e^"-e'3-«e-'"5,eflectcd. (3.10) 

There is no need to determine the phase factor e**", since 
it has no effect on the reflection probability |rp. 

A similar reasoning applies at the boundary x — L, 
where the matching condition reads 

= {vh/vd)^^^M'£, &tx = L. (3.11) 

Flux conservation requires that M' = '^^ e'" ''^ , 

with real parameters a',/3',7'. The value of a' is again 
irrelevant close to the Dirac point, because the spinor of 
the transmitted wave 

i^transmitted = tE^ ^^^^ (3-12) 

(with t the transmission coefficient) is an eigenvector of 
Ox- So 

^ ^transmitted — ^ ^6 (^^transmitted 5 (3.13) 

with a phase factor e*" that has no effect on the trans- 
mission probability 
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IV. TRANSMISSION PROBABILITY 

We consider the case W ^ L of a. wide and short crys- 
tal, when we may use periodic boundary conditions at 
y = 0,W for the Bloch waves \1/ oc e*^'^ ''. The transverse 
wave vector 5ky is then discretized at 6ky — 2im/W = 
Qn, with mode index n — 0, ±1, ±2, ±3, . . .. We seek the 
transmission amphtude i„ of the n-th mode. 

We first determine the transfer matrix Mn{x, 0) of the 
71-th mode $„(x)e**"^ through the photonic crystal, de- 
fined by 

$„(a;) = A/„(x,0)$„(0). (4.1) 

From the Dirac equation (12. 4p we obtain the differential 
equation 

—Mn{x,0)= ( — a, +g„(7, )M„(x,0), (4.2) 
ax \vd J 

with solution 

, ^ , 7 sinfc„a; / iSuj \ 
M„(x,0) = cosfc„x -I ax + qn<Jz ■ (4.3) 

kn \VD J 

We have defined the longitudinal wave vector 

fc„ = ^/iSco/vor - ql (4.4) 

The total transfer matrix through the photonic crys- 
tal, including the contributions l\3.6\i and (|3.1ip from the 



interfaces at a; = and x = L, is 

M = M''Hl„{L,0)M. (4.5) 
It determines the transmission amplitude by 

^0-r:)-C:)-(;;^)^^'(;:) 

2 2 

^ r = lT.T.^*r (4-6) 

where we have used the current conservation relation 

The general expression for the transmission probability 
Tn ~ l^riP is rather lengthy, but it simplifies in the case 
that the two interfaces at a; = and x ~ L are related by 
a reflection symmetry. For a photonic crystal that has an 
axis of symmetry at a; = L/2 both $(x) and ay^{L — x) 
are solutions at the same frequency. This implies for the 
transfer matrix the symmetry relation 

(JyMCTy = M^^ (TyM'cFy = M 

^ 13' = /3, y - -7, (4.7) 

and we obtain 
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^" cosh 2(3 — cos fc„L sinh 2(3 sinh 27 — ^' — - — sinh 2/3 cosh 27 



vnkn 



cos fc„Lcosh27 



qn sin knL , 



sinh 27 
I 



(4.8) 



For an ideal interface (when /? = = 7) we recover the 
transmission probability of Ref. Q . 

At the Dirac point, where Slo = Q ^ kn = iqn, Eq. 
reduces further to 



cosh^(q„L + 27) + sinh^ 2(3 sinh^(q„L + 27). (4.9) 



More generally, for two arbitrary interfaces, the trans- 
mission probability at the Dirac point takes the form 

^ = cosh2(/3 - (3') cosh^ C„ + sinh^(/3 4- (3') sinh^ fn, 



£,n = qnL + 7-7'- 



(4.10) 



While the individual T„'s depend on 7 and 7', this de- 
pendence drops out in the total transmission T„. 



V. PHOTON CURRENT 

The transmission probabilities determine the time av- 
eraged photon current I at frequency lojj + Slo through 
the photonic crystal. 



(5.1) 



n— — 00 



where /q is the incident photon current per mode. The 
sum over n is effectively cut off at \n\ ~ W/L ^ 1, 
because of the exponential decay of the T„'s for larger 
\n\. This large number of transverse modes excited in 
the photonic crystal close to the Dirac point corresponds 
in free space to a narrow range i5(/)Cia/L<Clof angles 
of incidence. We may therefore assume that the incident 
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FIG. 4: Same as Fig. [S] for /3 and 7 both nonzero. 
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FIG. 3: Frequency dependence of the transmitted current, 
for interface parameters f3' — l3, y' — —7. In the top panel 
we take 7 = and vary /3, while in the bottom panel we take 
/3 — and vary 7. The solid curves (/3 = 7 = 0) correspond to 
maximal coupling of the photonic crystal to free space. The 
curves are calculated from Eqs. (|4.8p and (jS.lf) . in the regime 
W/L ^ 1 where the sum over modes may be replaced by an 
integration over transverse wave vectors. 



radiation is isotropic over this range of angles S(j), so that 
the incident current per mode Iq does not depend on n. 

Since W/L ^ 1 the sum over modes may be re- 
placed by an integration over wave vectors, X]n=-oo ~^ 
{W/2'k) dqn- The resulting frequency dependence of 
the photon current around the Dirac frequency is plot- 
ted in Figs. [3] and 21 for several values of the interface 
parameters. As we will now discuss, the scaling with the 
separation L of the interfaces is fundamentally different 
close to the Dirac point than it is away from the Dirac 
point. 

Substitution of Eq. (|4.10p into Eq. (|5.ip gives the pho- 
ton current at the Dirac point, 

W 

I{Su; = 0)=IoTo — , 

_ arctan[sinh(/3 + P')/ cosh(/3 - f3')] 
°^ 7rsinh(/3-H/?')cosh(/3-/3') ' ^ ' 

independent of the parameters 7,7'. For two ideal inter- 



faces we reach the limit 



lim I(Suj 



0)//o 



1 w 



(5.3) 



in agreement with Refs. @, Eq. (|5.2p shows that, re- 
gardless of the transparency of the interfaces at x = 
and X = L, the photon current at the Dirac point is in- 
versely proportional to the separation L of the interfaces 
(as long as a <C i ^ W). 

As seen in Figs. El and [H the photon current at the 
Dirac point has an extremum (minimum or maximum) 
when either 7 or /3 are equal to zero. If the interface 
parameters /3, 7 are both nonzero, then the extremum is 
displaced from the Dirac point by a frequency shift Sojc- 
The photon current I{5lOc) at the extremum remains in- 
versely proportional to L as in Eq. (|5.2p , with a different 
proportionality constant Fq (which now depends on both 
fi and 7). 

The 1 /L-scaling of the photon current applies to a fre- 
quency interval \5uj\ "^vd/L around the Dirac frequency 
lod. For \5uj\ ^ wd/L the photon current approaches the 
i-independent value 



^00 ~ Ivi^- 



.W6uj 
ttvd ' 



(5.4) 



with rapid oscillations around this limiting value. The 
effective interface transmittance F is a rather complicated 
function of the interface parameters /3, 7, 7'. It is still 
somewhat smaller than unity even for maximal coupling 
of the photonic crystal to free space (F = 7r/4 for = 
7 = 0). 



VI. CONCLUSION 

While several experiments [13, El have studied two- 
dimensional photonic crystals with a honeycomb or tri- 
angular lattice, the emphasis has been on the frequency 
range where the band structure has a true gap, rather 
than on frequencies near the Dirac point. Recent ex- 
periments on electronic conduction near the Dirac point 
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of graphene have shown that this singularity in the band 
structure offers a quahtatively new transport regime . 
Here we have explored the simplest optical analogue, the 
pseudo-diffusive transmission extremum near the Dirac 
point of a photonic crystal. We believe that photonic 
crystals offer a particularly clean and controlled way to 
test this prediction experimentally. The experimental 
test in the electronic case is severely hindered by the 
difhculty to maintain a homogeneous electron density 
throughout the system [l^. No such difficulty exists in 
a photonic crystal. 

If this experimental test is successful, there are other 
unusual effects at the Dirac point waiting to be observed. 



For example, disorder has been predicted to increase — 
rather than decrease — the transmission at the Dirac 
point [3, [iM [lB|- Photonic crystals could provide an 
ideal testing ground for these theories. 
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